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Abstract 

We review some recent theoretical progress in the understanding of weak radiative 
B meson decay, and discuss the implications of present and improved measurements 
of b — ► S7 for supersymmetric models. 
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1 Introduction 



Measurements of flavor changing neutral current processes provide some of the most stringent 
experimental tests of the standard model and its extensions. In the standard model, the 
smallness of quantities such as the K — K mixing amplitude is linked to the fact that all 
such flavor changing amplitudes vanish at tree level, and can proceed only at higher orders in 
the weak interactions. Because the standard model predictions are so small, these processes 
provide an excellent testing ground for extended models possessing new sources of flavor 
violation. 

Unfortunately, it is not always a simple task to convert the underlying quark-level weak 
amplitudes into robust predictions for the observed hadronic quantities. Strong interaction 
effects can significantly modify the quark model estimates of decay rates. Luckily, in the case 
of the B mesons, the large mass m b of the b quark makes it possible to systematically compute 
the effects of strong interactions in a double expansion in a s {m b ) ~ 0.2 and AQ C D/ m b, where 
Aqcd ~ 200 MeV is an energy scale typical of hadronic quantities such as the hyperfine 
splitting in the heavy mesons. Hence the rare B decays B — > X s ^j and B — > X s £ + £~ with 
£ — e, /i,t are especially suitable as tools for testing the standard model and its extensions. 

Over the last year or so, there has been some theoretical progress pertaining to the 
inclusive b — > S7 transition. With respect to the perturbative calculations, the complete 
next-to-leading order (NLO) QCD corrections have been completed, resulting in a significant 
reduction of the theoretical error. In addition, an interesting new class of non-perturbative 
corrections to the rate have been identified. These effects are of order AQ CD /ml, with m c 
the mass of the charm quark. Such effects are potentially important, of order 10% or more. 
Furthermore, they appear, on the face of it, to violate the expectations of heavy quark 
effective theory, which predicts that all such corrections are suppressed by inverse powers of 
the b quark mass. 

In the following section, we review the standard model prediction for the branching ratio 
B[B — > X s j). In Sec. 3, we discuss the impact of supersymmetry on the rate, and estimate 
the range of superpartner masses that might be probed by improved measurements. Sec. 4 
concludes. 



2 Standard Model Theory 

2.1 Status of Perturbative QCD calculations 

In the standard model, the b — > 57 decay is mediated by the W exchange diagrams shown 
in Fig. 1. Ignoring QCD corrections, these diagrams result in an effective Lagrangian of the 
form 

C^ S1 = ^0V*V tb C 7 s L a^b R F, u , (1) 

where C 7 contains the dependence on the W and top quark masses, and is given by 

3x^ 2,x^ 8x^ -\- 5x 7x o/-. r o z„\ 

C ^K^W^ X - 24(,-l)3 ■ * = ™?/«S, (2) 
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Figure 1: Feynman diagrams for b — ► 57. 



Perturbative QCD corrections can be included by use of the renormalization group equations, 
which enable us to scale the coefficient C-j from a high scale fiw ~ Mw to a low scale //j ~ m^. 
In this way, the large logarithms of order a s log(M^ / m&) can be resummed into the coefficient 
function CV(/x), which now depends on the renormalization scale fi as well as the top and 
W masses. The leading order QCD corrections have been known for some time 0, and 
the NLO corrections have recently been completed 0. In calculating the branching ratio, 
it is convenient to normalize to the semileptonic width. This eliminates a large uncertainty 
arising from the usual factor mf appearing in weak decay rates. The ratio is given by 

B{B -> X c ev e ) \V cb \ z ttJ-(z) 

Here z = ^c, P oic/ m fe,poie5 an d /( z ) i s the phase space factor for the semileptonic width. F 
includes the QCD correction to the semileptonic width, while D is the QCD corrected value 
of CV(/x), and A gives the contribution from certain infrared logarithms appearing in the 
b — > S7# rate. The full expressions for /, F, D, and A can be found in Ref. [0. Finally, the 
quantity <5np is the the correction due to non-perturbative binding effects, which we shall 
return to below. The important point here is the dependence on the renormalization scale 
fib- Physical quantities should not depend on fib, although at any finite order of perturbation 
theory, some residual dependence on fi b remains. Typically /x& should be chosen of order m b , 
but the variation of the prediction for a range of //& ~ m b /2,2m b indicates the uncertainty 
resulting from higher order corrections. In leading order, this uncertainty is roughly 30%; 
the inclusion of the NLO corrections reduces this to about 4%. 



2.2 Non-Perturbative Corrections 

Apart from the perturbative corrections discussed so far, the rate is also modified by non- 
perturbative binding effects. These effects can be systematically calculated by means of the 
heavy quark effective theory || together with the operator product expansion (OPE). Since 
the energy release in B decays is large compared to the typical strong interaction scale Aq CD , 
the corrections can be organized as a series in powers of AqcD/^fe- The leading corrections 
are governed by the quark 'pole' mass m b pole , a quantity A that measures the mass of the 
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Figure 2: Triangle diagram for the b — > s^g process. We have omitted a second diagram 
with the gluon and photon interchanged. 



light degrees of freedom in the meson, and finally two quantities of dimension mass 2 that 
measure the Fermi motion of the b quark and the hyperfine splitting between the B and B* 
mesons, respectively. Ai and A2 can be expressed in terms of B meson matrix elements as 

(B\b(tD) 2 b\B) (B\b(ga^G, u )b\B) 
Al " 2M~ B ' As " UM~ B • (4) 

The quantities Ai,2 ~ Aqcd and A ~ Aqcd suffice to compute the non-perturbative effects 
through order 1/m 2 . For the b — > 57 decay, the rate is given by @] 



r(a-x. 7 ) = r (i + \^), (5) 



where r is the quark- level decay rate. Including the 1/m 2 contribution to the semileptonic 
decay rate ||, we find that the ratio R in Eq. @ is enhanced by roughly 1% through non- 
perturbative effects. Note in particular that when the width is expressed in terms of the 
quark pole mass, the leading corrections are at order 1/m 2 and hence are expected to be 
quite small. 

It has recently been shown that the inclusive B —>■ X s 7 rate is subject to corrections at 
order 1/m 2 |5[]. On naive dimensional grounds, these corrections could be sizable. Indeed, 
the scale of these effects is set by quantities like A2/m 2 ~ 0.08, so shifts of 10% or more 
are possible. The 0{l/m 2 c ) contribution to the inclusive rate comes from an amplitude 
involving the b — > ccs transition, where the cc pair annihilates into the final state photon 
after absorbing a soft gluonQ The resulting loop diagram is shown in Fig. 2. Similar loops 
involving u and t quarks can be safely neglected. The dependence of this diagram on the 
charm quark mass can be deduced using dimensional analysis || . In terms of the photon 
and gluon momenta fc 7 and k g , the diagram depends on the invariants /c 2 , fc 2 , fc 7 • k g , and 
m c . In the limit of an on shell photon and a very soft gluon, all of the invariants involving 
four momenta vanish, and we conclude that the amplitude depends on m c only. Explicit 



2 Similar corrections of this type involving hard gluons have been considered in Refs. Q 
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evaluation yields a new contribution to the weak effective Lagrangian which mediates the 
b — > S'yg transition: 



eQ c C 2 



X" 



V2G F V; s V cb (s L r^b L )^-G: x d x Ft 



3m;? 



(6) 



where Ci ~ 1.02 [[| is the coefficient of the b —>■ ccs four- Fermi operator. Including this term 
in the OPE for the inclusive B — ► X s 7 rate, we find a correction of order 1/m 2 H. Efl: 



r(5^X s7 ) " 9C 7 m2 



+0.025. 



(7) 



In fact, the approximation leading to Eq. ([/]) is not as robust as one might hope 10]. The 
evaluation of the triangle diagram involves an integral of the form 



I(ky, kg) 



24 / xdx 



l-x 



ydy- 



1 



m 2 — k 2 x(l — x) — k%y(l — y) 



(8) 



< mi 



which reduces to 1/m^ in the limit k£, k 2 , ■ k g <C m\. In actuality, we have k 2 

Aq CD <C m 2 c \ however k y ■ k g ~ mfcAQCD ~ Hence the result ([FD, which 



and k 2 g 



involves setting /c 7 • k g to zero, is potentially subject to large corrections. These corrections 
manifest themselves as a series of higher dimension operators in C b ^ s ~t9. Operators of di- 
mension (8 + n) in this series will contribute to the inclusive rate by an amount proportional 
to (mbAQ CD / 'm 2 c ) n '. Since (mbAq CD / 'm 2 c ) ~ 0(1), these subleading contributions are not 
suppressed by any small number. 

The first few higher dimension operators in the b — > s^g Lagrangian can be found by 
expanding the charm loop integral in powers of fc 7 • k g /m 2 c . The effect of these terms on 
inclusive rate can be found as before using the OPE: we have [ PU[ | 



6T(B -> X sl ) 



Co 



9C 7 12M B 



[B\gb 



—a uX G + ——a va [%D )G 



mi 



15 mi 



3 m\ 
700^! 



a vX {iD) 2 G' x + o va {iD x ,iD a }G 



b\B). 



(9) 



Here we have specialized to the rest frame of the B meson, where to leading order Pf, = 
(m&, 0). The first term in square brackets gives the leading contribution of Eq. (0), while 
the subsequent terms give the first two subleading corrections. In fact, the first subleading 
term vanishes at this level of approximation by virtue of the QCD analogue of the Maxwell 
equation V x E = — B, together with the heavy quark equations of motion [ID|. The next 
subleading term is certainly negligible, owing to its small 3/700 coefficient, which arises from 
solid angle factors in the photon phase space integration as well as small coefficients in the 
Taylor expansion of the loop integral. The smallness of these corrections leads us to believe 
that the leading ~ 3% result of Eq. (^) is in fact quantitatively correct. 

It is interesting to ask whether non-perturbative 0(l/m 2 c ) effects of this type might be 
important in other B decays |IT| . The impact of the charm loop onB-* X s e + e~ has been 
estimated || Yl\, and found to be similar to that in B — > X^, of order a few percent. It is 
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unlikely that corrections of this type could have a significant impact on the B semileptonic 
branching ratio or the ratio of the B° and A& lifetimes. A comparison of A& — > X s , -f with 
B — ► X s 7 is one way to test for spin dependent effects of this type |J . Since the light degrees 
of freedom in the B carry spin 1/2, while those in the A& carry spin 0, we expect that the 
charm loop will have no effect on A& decays. 

Finally, we should point out that the 1/m 2 contributions are not in conflict with heavy 
quark effective theory. Indeed, in the HQET limit m& — * oo with m c fixed, the corrections 
are not proportional to l/rn| and have a rather different form |T(J. The presence of these 



effects can be traced to a breakdown of the OPE in processes where the final state photon is 
radiated from light quarks, rather than being produced directly in the b — > s transition 11. 

2.3 Standard Model Prediction 

Assembling the ingredients of the previous subsections, the standard model prediction for 



the branching ratio is 0, 113 



B(B -> X sl ) = 3.48 ± 0.13 ± 0.28 x 10" 4 , (10) 

where the first error is due to renormalization scale dependence and the second is due to 
errors on the input parameters, which have been added in quadrature. A detailed accounting 
of the errors can be found in 0, [TB|]; the dominant errors at present are due to the ratio 
of the b and c quark pole masses (5.3%), the semileptonic branching ratio (4%), and the 
renormalization scale dependence (4%), with smaller ( ~ 3%) errors from a s (Mz), CKM 
angles, and the top quark mass. One promising way to reduce the theoretical error is through 
better determination of the quark pole masses. Analyses of the lepton energy spectrum or 
the hadronic invariant mass spectrum in semileptonic decays have been proposed as means 



of extracting the pole masses |14j . In addition, measurements of the photon energy spectrum 
in B — > Xsj can be used to provide largely independent constraints on the HQET input 
parameters including the pole masses. 



The prediction (|1(]) is to be compared with measurements by CLEO [|T^] and ALEPH 



17 : 



B{B^X sl ) = 2.32 ± 0.57 ± 0.35 x 10~ 4 (CLEO) (11) 
= 3.38 ± 0.74 ± 0.85 x 10" 4 (ALEPH). (12) 



For reasons discussed in [ 13 1 , it is important to rely as little as possible on theoretical models 
of the photon energy spectrum in measuring the branching ratio. Indeed, measurements of 
the spectrum can be used to refine our theoretical understanding of the B system. As 
an example, the photon energy spectrum can be related to the lepton energy spectrum in 



inclusive semileptonic b — ► u decays [18 



3 b — > 57 and Supersymmetry 

The impact of new physics, and in particular supersymmetry, on flavor changing neutral 
currents is the subject of a vast and growing literature. The case of b — > S7 has been studied 
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extensively f[5 |. In this section we review some of the salient features of supersymmetry as 
it relates to flavor changing neutral currents in general and the b — > 57 decay in particular. 
Supersymmetry (SUSY) predicts the existence of a scalar partner for each standard model 



fermion, and a fermionic partner for each gauge or Higgs boson [20]. The left- and right- 
handed fermions have separate scalar partners, which can have different masses and can mix 
with one another. In addition, softly broken supersymmetry requires at least one additional 
Higgs doublet in order to generate masses for both up- and down-type quarks. The Higgs 
spectrum then consists of two neutral scalars, a neutral pseudoscalar, and a charged Higgs. 
When the electroweak symmetry is broken, the charged gauginos and Higgsinos mix to form 
two Dirac 'charginos,' while the neutral gauginos and Higgsinos mix to form four Majorana 
'neutralinos.' 

Flavor changing neutral currents are problematic for general supersymmetric models. 
This is a consequence of flavor mixing among the squarks. If the squarks are non-degenerate 
and strongly mixed, the couplings of the gauginos and Higgsinos will not conserve flavor. 
Loops involving (for instance) squarks and gluinos can then give large contributions to flavor 
changing processes. Hence there is a potential conflict between weak-scale supersymmetry 
and small flavor changing neutral currents |2l] . 



These problems are alleviated if it is assumed that the squarks are degenerate in mass. 
This is the case, for instance, in supergravity |22| and gauge-mediated [^3] models of super- 
symmetry breaking. In this case, the gluinos and neutralinos have flavor-diagonal couplings, 
while the chargino couplings reduce to the usual standard model CKM matrix. Renormal- 
ization effects break the degeneracy of the squarks. Typically, the squarks will be degenerate 
at a high scale (~ 10's of TeV to Mpi anc k, depending on the model) where supersymmetry 
breaking is communicated to the standard model particles. When the squark masses are 
renormalized down to the weak scale, the squarks of the third generation are lighter than 
the others. Even in the presence of this mass splitting, the gluino, neutralino, and chargino 
couplings are not strongly modified in comparison to the case of exact degeneracy. Tests of 
'minimal' models of gauge mediated supersymmetry breaking based on b — > 57 have been 



discussed in [24 



To illustrate the sensitivity of b — > 57 to high mass superpartners, we will consider 



a relatively benign scenario with the following parameters [25]. All squarks, except for the 



lighter top squark t 1; have a common mass of 500 GeV. The light stop mass is set to 200 GeV, 
and is taken to be a maximal mixture of and Ir. The ratio of the vacuum expectation 
values of the two Higgs doublets, tan @, is equal to 3, and the charged Higgs mass is 600 GeV. 
Finally, we will assume that the neutralinos and gluinos have flavor diagonal couplings, while 
the chargino couplings are proportional to the standard model CKM matrix. In this limit, 
the dominant contributions to b — > 57 arise from loops involving W and the top quark, the 
charged Higgs and the top quark, and the charginos and the top squarks. The charged Higgs 
enhances the rate, while the charginos can enhance or suppress the rate. Consequently it is 
not possible to derive model-independent bounds on superparticle masses, although in non- 
supersymmetric two-doublet models, the lower bound on the charged Higgs mass turns out to 



be around 340 GeV pq] . To compute the rate, we must specify two final parameters: /i, the 
Higgsino mass parameter, and M 2 , the SU(2) gaugino mass. These parameters determine 
masses and couplings of the charginos. Fig. 3 qualitatively illustrates the dependence of 
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Figure 3: Contours of constant chargino mass (solid curves) and constant B(B — > X s ~f) 
(broken curves). Both axes and the chargino masses are in GeV units. See text for details. 



B(B — > X s j) on fi and M2 for this choice of parameters^. For clarity of presentation, we 
have shown only fi > 0; with this choice of stop mixing, the charginos enhance the rate for 
fi < 0. The solid curves are contours of constant chargino mass, along which the lightest 
chargino has mass 100, 200, 300, and 400 GeV. The broken curves are contours along which 
the ratio r = B(b — ► s7)susy/-E>(& — * s7)sm * s constant, with values 0.3, 0.5, 0.7, and 0.9. 
We see that deviations from the standard model at the 10-20% level are possible for chargino 
masses in the 300-400 GeV range. Consequently, improved measurements will continue to 
narrow the allowed region of SUSY parameter space. 



4 Conclusions 

Theoretical understanding of B — > continues to improve. The completion of the NLO 
QCD corrections has dramatically reduced the uncertainty on the standard model prediction. 
The non-perturbative corrections of order 1/mj* are well understood, and the corrections of 

3 We caution the reader that while Fig. 3 includes all NLO corrections for the standard model case, the 
SUSY prediction does not include NLO matching conditions for the Wilson coefficients. These have not yet 
been evaluated. This approximation is adequate for our qualitative purposes. 
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order while unexpected, are both small and calculable. The sensitivity of B — > to 

new physics is quite remarkable, and we have reason to hope that improved measurements 
will continue to add to our knowledge of short-distance physics. Improved measurements of 
the photon energy spectrum will also be an aid in our understanding of hadronic physics. 
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